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^ ' Abstract 



(N 



A detailed study of the heavy quark h = c,b, . . . contributions to deeply virtual Compton 
scattering is performed at both the amplitude and the cross section level, and their phe- 
nomenological relevance is discussed. For this purpose I calculate the lowest order off- forward 
] photon-gluon scattering amplitude with a massive quark loop and the corresponding hard 

' scattering coefficients. In a first numerical analysis these fixed order perturbation theory 

r~| ' results are compared with the conventional intrinsic "massless" parton approach consider- 

J^, ing generalized parton distributions for the heavy quarks. The differences between these 

Q^l two prescriptions can be quite significant, especially at small skewedness where the massless 

^ ■ approach largely overestimates the deeply virtual Compton scattering cross section. 



■ 1 Introduction 



In deeply virtual Compton scattering (DVCS) El 0] a highly virtual photon, generally radi- 
ated from a charged high-energy lepton, converts to a real photon by scattering on a nucleon 
target that remains intact. The phenomenology of DVCS looks very promising with first data in 
different kinematic regions being available from several groups at fixed target jH ^ and collider 
experiments [HI E]. On the theoretical side the amount of uncertainties in the predictions for 
DVCS could be reduced by incorporating next-to-leading order corrections in perturbative quan- 
tum chromodynamics (QCD) |Hl El HH] and power corrections to the leading twist-two results 
starting at the twist-three level [TTl |T2l UH] • 

The main motivation for studying DVCS is the access to new non-perturbative information on 
the structure of the nucleon carried by generalized parton distributions (GPDs) [TU[TJ|2l [13 Ellin], 
which will deliver insight into the spin structure of the nucleon, in particular the quark and gluon 
spin and orbital angular momentum contributions, PP and into the nucleon structure in three 
dimensions ^Zj. An excellent survey of the theory of GPDs and their experimental accessibility 
can be found in a recent comprehensive report jl8j . 

It is reasonable to inspect the inclusive proton structure function F2 in deep inelastic scattering 
(DIS) due to its close connection to the general Compton amplitude via the optical theorem. In 
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piP.S) p{P',S') 
Figure 1: General kinematics of deeply virtual Compton scattering on a proton. 

the kinematic region relevant for DVCS at the colhder experiments the charm contribution F2 
to the proton structure function has been found to be up to 20% O 120] • If one presumes this 
same ratio for the corresponding form factor of the DVCS amphtude this results in a charm 
contribution of about one third to the DVCS cross section (amplitude squared). Hence a proper 
treatment of charm in DVCS can be quite important. 

Generally, heavy quarks h = c,b,t are characterized by ml 3> ^qcd whereas light ( "massless" ) 
quarks have mf^^^ ,, <^ ^qcd- fixed order (FO) perturbation theory, heavy quarks contribute 
in lowest order to the plioton-gluon scattering (PCS) Y9 ~^ 19 DVCS subprocess where they 
appear as massive internal fermion lines. As is well known, such a fixed order calculation correctly 
describes the threshold region s > Ami in DIS [11112011^11221, and its perturbative stability has 
been demonstrated in Ref. [22] even up to s ~ 10^ GeV^. In the conventional parton model 
approach potentially large logarithms log"(m^/s) appearing in the fixed order PCS result are 
resummed to all orders and absorbed into intrinsic "massless" parton (MP) distributions for the 
heavy fiavors j24j . 

In the next section, after giving a rough overview of DVCS in the leading twist approximation, 
the technicalities of incorporating heavy quark contributions in the DVCS amplitude are described 
and the need for massive hard scattering coefficient functions is explained. In Sec. [21 a systematic 
way to calculate the lowest-order of PCS with a massive quark loop is introduced whose relevant 
results are presented in Sec. 01 This completes the requisites for a first numerical analysis of 
heavy quark effects in DVCS where the numerical results at the amplitude as well as at the cross 
section level are discussed in Sec. El Finally, in Sec. [HI the conclusions are drawn and an outlook 
is given. 

2 Leading twist DVCS amplitude 

First, some general properties of DVCS have to be reiterated where mainly the notations of 
Refs. [23 E] are used. The general kinematics of DVCS on a proton are shown in Fig. ^ with the 
masses defined to he := — and := = P'^. An independent set of momenta is given 
by 

P:=(P + P')/2, g:=(g + g')/2, A:=P'-P (1) 
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that leads to the following five Lorentz invariants 

g^:=-g^ x:=^, ^= ' ^ - ' := ^ = - t/4 (2) 

2F-q 2F-q 

where for the moment the outgoing photon is also considered virtual. Only the first three in- 
variants survive the generalized Bjorken limit Q"^ ^ oo with x, ^ fixed in which the first two are 
straightforward generalizations of corresponding invariants in DIS. The additional scaling variable 
^, the so-called skewedness, is a measure of the difference between the two photon virtualities. 
Furthermore, the following invariants prove to be helpful for the representation of the analytical 
results in Sec. HI 

s := {P + qf = g^(-l + 1/x) + W , u:={P-qf = g2(_i _ ij^) + ^2 ^ 
q^ = {q + A/2f = g^(-l - e/x) + t/4 , q'' = (q - Ajlf = g^(-l + ^x) + t/4 (3) 

where and t/4 are neglected in the generalized Bjorken limit. In DVCS with a real photon 
in the final state only one scaling variable remains 

x^i = x{\- t/AQ^) , = QV2 + t/4 ^ , (4) 

which is chosen to be ^ to avoid confusion with the usual Bjorken variable xb ~ 2^/(1 + C.) in 
DIS, and it is more appropriate to use Q^. 

The dynamics is contained in the DVCS amplitude 

-iT^^ = j e'^-'{P'S'\T[r{z/2)J''{~z/2)]\PS)d^z (5) 

that is given by the time-ordered product of two currents which are related to the electromagnetic 
current by J^^iz) = eJ^{z). The leading twist contributions, apart from those that arise due to 
photon helicity fiip which are beyond the scope of this work, are contained in the following two, 
transversal photon spin conserving, form factors 

= E ^A^r + E ^^A^r t^+--- (6) 

A=± A=± 

whose exact Lorentz structures will be defined when they are explicitly needed in the next section. 
For now it is sufficient to note that they are equivalent to the usual ones in the generalized Bjorken 
limit, e. g. the tensors i^u {Eq. (9) in Ref.jH]}. These two form factors can be further decomposed 
with respect to their Dirac structure 



Tt = U{P', S') 
Tt = U{P', S') 



2p.q ' AMP-q _ 

q-A 



2P-q 4MP-q 



U{P. S) , (7) 
75f/(P, S) (8) 



in analogy to the form factors of the electromagnetic current between two nucleon states. 
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On the basis of the factorization theorems proven in Refs. fl^ |23 I2Z| each of these form 
factors can be expressed as convolutions of GPDs with perturbatively calculable hard scattering 
coefficients such as 



a=q,g-^-^ ^ b=q' 



(9) 



Equivalent formulas can be obtained after the replacements Trp^,Trp2 and H^,E^, 

trp^^. The GPDs in Eq. © are defined in 



or Trr 



^T,A5 ^T,P5 fa ~^ fa}Ha,E^, and trp . ^ 



such a way that they have the closest possible connection to the unpolarized and polarized DIS 
parton distribution functions, that is Hg g{x, 0, 0, /i) = q, g{x, /i) and H^ g^x, 0, 0, /i) = Sq, 5g{x, fi), 
respectively. The factorization scale fi is as usual assumed to be equal to the renormalization 
scale with the common choice fi^ = Q"^. The hard scattering coefficients describe the scattering 
process 7*0 ^ 7a of the intrinsic parton a with the coupling to the photon being mediated by 
the quark b and accordingly define the sums in the factorization theorems. 

Crossing symmetry of the Compton amplitude allows to rewrite the hard scattering coefficients 
in the following form 



(10) 



where the definition of the coefficient function (C^p is not unique since odd (even) pow- 

ers of z drop out in the relevant combinations. Eventually, the perturbative expansion of the 
coefficient functions is written as 



'11^ 



The MS coefficient functions with exclusively massless parton lines are available up to next- 
to- leading order [23 I2H] • For definiteness and since some of them are needed for comparison in 
Sec. m they are explicitly given by 



C^^Uz,Q\fi) = C'^nz,Q'\f,) + 6,a 



a 



(0) 



-6. 



ba 



9 



z- 1 
1 



Z{1 

3 

z + 1 



It 



log 



2z 



2{z+l] 



log 



z-l 
z-l 



4 



2z 



(1) 



3 + 21og 

6 — loj 

z-l 



z-l 
z-l 



log^ 



2z 



2z 



2z 



- - 2 log — 



z+l 



log 



2z 



log 



z - 1 



2(z+l)2 ° 2z 



4 , z 

- 2 + log - 



z-l 



2z 



(12) 



(13) 
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None of the coefficient functions that contain internal massive quark hues have been available so 
far and it is the aim of the next two sections to obtain these in lowest order. 

In the general form of the factorization theorem in Eq. Q the inclusion of heavy quarks is 
simply done by appropriate specification of the GPDs and the accompanying hard scattering 
coefficients. In the fixed order perturbation theory approach one has GPDs for the three light 
flavors M, (i, s and for gluons. In the hard scattering coefficients the light flavors are considered 
massless whereas the masses of the heavy quarks are kept. Explicitly, in leading order one has 




h=c,b,t 

with an analogous formula for the helicity dependent form factor Trp. The factorization scale of 
the fixed order part can be chosen independently and is preferably /i/j = 2mh |2S]- Apparently 
there seems to be a mismatch in the order of the hard scattering coefficients as the lowest order 
massive ones start at 0{as) that can be resolved by comparing the respective magnitudes which 
turn out to be partly of the same size. Additionally this is in line with the "massless" parton 
approach where heavy quarks contribute on the same level as the light quarks. In this latter case 
one considers GPDs for all quark flavors where the intrinsic heavy quark GPDs are generated by 
the usual massless evolution equations [HJ El 121] with the boundary conditions, in analogy to 
Ref. El, 

fhi^, tt,fi <mh) = , fh{x, ^, t, /i < ruh) = (15) 

and the massless coefficient functions are used throughout. Specifically in Eq. (jl4p the ffist sum 
is extended to all flavors and the second sum with the gluon GPD is absent. 



(14) 



3 Photon-gluon scattering 

The lowest-order Feynman diagrams that contribute to photon-gluon scattering are shown in 
Fig. 121 The diagrams with the quark direction reversed are not shown separately but are indi- 
cated as antiquark loops. Even though the reversion leads to topologically different diagrams in 
configuration space, their evaluation leads to the same result and is taken into account by an 
overall factor of two in the following PGS amplitude in D dimensions 



t 



X 2Tr 



+ 



+ 



(27r)^ 



+ f + 4' + mh)j''{it + ^ + mh) 



Dk{Q, nih) Dk{p', rrih) Dk{p' + g', rrih) Dk{p, rrih) 
7°(^ + mfe)7^(^ + / + mfe)7^(^ + j) - j' + mh)Yi^ + i> + ^h) 

Dk{0, rrih) Dk{p', mt) Dk{p - q', rrih) Dk{p, rrit) 
7"(^ + mh)Y{^ + + mh)-f^{lt + f + ^' + mh)-i^{f + + m^) 

Dk{0, rrih) Dk{q', rrih) Dkip' + g', rrih) Dk{p, rrih) 



(16) 
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9aip,£") 9b{p -.e'^*) gaip,£") 9b{p',£'^*) gaip,£") 9bip',e'!^*) 



(a) (b) (c) 

Figure 2: Lowest order Feynman diagrams for photon-gluon scattering with color indices a,b,i. The 
diagrams that have the heavy quark hi direction reversed are indicated by h^. 



with Dk{p, rufi) := {k — ml + ie. The calculation of this amplitude was carried out with the 
help of Mathematica [30] that provided an adequate handling of the excessive algebra and the 
Tracer jST] package that was used to evaluate the trace and to contract Lorentz indices. The 
individual steps that lead to the analytical results are detailed below. 

The integral over the loop momentum can be systematically reduced to scalar integrals 

J d^k[DkiO,mo) Dk{pi,mi) ■ ■ ■ DkipN,mN)]~^ (17) 

by a Passarino-Veltman decomposition [32] which is described in a form more suitable for imple- 
mentation in computer algebra systems in Ref. j^SI- The whole PGS amplitude in Eq. is 
finite whereas the separate scalar integrals can be divergent. Therefore the general D dimensions 
have to be kept to regularize any divergences during intermediate steps. Explicit formulas for the 
scalar integrals are given in Refs. [3U EHl tHH] in which they are expressed in terms of complex 
(di)logarithms whose arguments are generically built up by 



with 13, — ^ 



I3s+1 ^ \ s + ie 

if all masses are equal as it is the case for the lowest order PGS amplitude. 

The immediate application of the Passarino-Veltman decomposition to Eq. ()16|) is in principle 
possible but not practical. The amount of work is considerably reduced if one gets rid of the four 
free Lorentz indices first. Inspired by Refs. (SHI IHZI covariant tensors with an obvious physical 
meaning are constructed that appear as independent Lorentz structures in the PGS amplitude 
and equally act as projectors onto them. The three covariant four-vectors 

.= r^^^ - i■r]7]^' - y'-gg'' ^ _^ _^ r/-gg^ - g'r^^ ^^^^ 

with r:=4^^, V' = l, 
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7]' defining the direction of 62, and definite parity serve as a basis. One easily proves that they 
reduce to the usual linear polarization vectors of the incoming photon in the frame ff = 0. Now 
the sought-for tensors are simple products of these polarization vectors and analogous ones for 
a, P, with the sensible choice r] ^p + q = p + q and rj' = q. Note that the totally antisymmetric 
epsilon tensor always appears in pairs as the PCS amplitude is parity-even. These products are 
easily generalized to D dimensions if they are expressed by metric tensors. Physically, this means 
a sum over polarization vectors that are perpendicular to the scattering plane. In principle this 
sum has to be normalized by a factor 1/{D — 3) which is, however, not necessary for the finite 
PGS amplitude in Eq. (fT^ . 

The complex circular polarization vectors used in the expansion of the Compton amplitude 
in Eq. (0) are given by 

e^:=T^{e1±zeii) , ^ := (20) 

From a parton model point of view it looks advantageous to use rj ^ p because in this case 
these polarization vectors are invariant under the replacement p —>■ yp. Anyway, both choices are 
equivalent on the leading twist level and lead to the same hard scattering coefficients that are 
presented in the next section. 

In order to extract the leading twist contribution, one has to consider the limit t —>■ where 
some care has to be taken since factors up to that appear in denominators which originate 
from the normalization of the polarization vectors have to be canceled thoroughly. In this limit 
the dilogarithms drop out completely. 



4 Analytical results 

At the leading twist-two level the general Compton amplitude with virtual photons in the initial 
and final state receives contributions of gluons inside the proton to the following form factors 



X=± A=± 

+ E ^-vr E ^^-vr ?t + • • • (21) 

A=± A=± 

where the three additional terms in comparison to Eq. (jH)) are related to longitudinally polarized 
photons (L) and to photon helicity fiip (T). Due to the real photon in the final state T-^ does not 
appear in the DVCS amplitude. From the calculation of PGS that has been sketched in the last 
section it is possible to obtain the lowest-order massive hard scattering coefficients appearing in 
factorization theorems for these form factors like in Eq. Q. Now it will be indicated how they 
are extracted from the amplitude in Eq. 

Generally, the squared heavy quark charge already appears in the factorization formula, the 
QCD coupling constant is extracted in accordance to Eq. and one gets an additional factor 
1/[2(1 — ^^)] from the normalization of the gluon GPDs |TH1 EH]- The massive hard scattering 
coefficients relevant for the numerical analysis in the next section are contained in the gluon 
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polarization vector sums ^'\^T ^T) resp. Ae^e"* for Trp and read 



+ 2 [(1 - 2xf - + 4f7a;(l - x)] fogrj 

- [(x/^ + 1 - 2f7x/^) (1 + - - 4f7x^ - 8ffx^] fog^ 

- 2 (x/e + l){l + Ax^-e- 4?7<) fogv +(x -x)} 
= 2iV^ey [ - 2x + fog^ r, +2 (3 - 4x + (3, fog r, 

+ 2 (x + (fog' r,2 +4 fog r,2) - {x ^ -x)] (22) 

with f/ := ml/Q'^. The corresponding results for massless quarks j2Sll2H] are obtained for f/ — 
with appropriate subtraction of the logarithmic divergences ~ log(m/i//i). It should be remarked 
that the limits have to be taken in this order, namely, first t — >■ and then — 0, to reproduce 
the massless results. 

Furthermore, the contraction of the gluonic Lorentz indices analogously to Trp gives the con- 
tribution of longitudinally polarized photons 

= ""(1 _y)2 [2f7x2 (logV,-logV,2) +2x(l -x)/5,logr, 

+ {x/^ + 2x^ + x^) Pg2\ogrg2+{x ^ -x)] . (23) 

The "massless" limit can be achieved in the same way as mentioned above but will not be explicitly 
stated since it is not needed at present. 

As a byproduct of the calculation of the PGS amplitude with the method presented in the 
previous section one also gets massive expressions for the hard scattering coefficients in factor- 
ization theorems with generalized helicity-fiip, or tensor gluon distributions, respectively, whose 
exact definitions and applications are given in Refs. [HHlEnilini- They have to be projected from 
the gluon polarization vector sums for and Xe'^^e'^* for T-^ and are given 

for completeness 

4m = = (^^^^2)2 [2^(1 + (log' - log' - 2^a;(l + log r, 

- {x + C) {x-^- 2f]x) {(3^ hgrs - 13^2 \0grq2) + {x -x)] 

The limit ^ is finite as expected because the generalized tensor gluon distributions do not 
mix with any quark distributions and reproduces the "massless" results of Refs. |3? 
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Two further limits provide additional checks for the mass dependencies in Eq. fl22|) . On the 
one hand the contributions of heavy quarks do vanish for infinitely large masses in accordance 
with the decoupling theorem |1T] . On the other hand the imaginary parts of Eqs. (j22|) and (j23|) 
in the forward scattering limit are given by 



im . > 



4r/) 



X { [x^ + (1 - xf + 4?7x(l - x) - Sr/V] log(- r^) + [(1 - 2xf + 4?7x(l - x)] /5 J 

47]) [2r7x^ log(— Ts) + x(l — x) /?^] 



-Im t« i=^4T^^(l/x-l 

TT ' 



— Im t. 



7(1) C-0 



TT 



Tp 0(l/x - 1 - 4r7) [(1 - 2x) log(- rg) + (3 - 4x) (3^ 



(25) 



with 7] := m\/Q'^. They coincide with the corresponding quantities of the DIS structure functions 
{F2 — Fl)/x, Fl/x jl2j, and 2gi [43j in view of the factorization theorem Eq. (jU)) as demanded 
by the optical theorem. 

Finally the relevant results for the numerical analysis in the following section are the massive 
DVCS coefficient functions as defined by Eq. (fTO| 



2(^+1)2 



z - 1 



[Ps \ogrs-(3g2 log r 52) 



+ 



1 + 8t]z 
l-2r] 



1 - (1 + 4:r])z 
1 + (3 + 877)^2" 



log^ rs 



log^ r 



2(z + 1)2 [ Z 



^ {pg log rs - Pg2 log rg2) + log^ rs - log^ rg2 



(26) 



that are achieved by enforcing the DVCS kinematics x —>■ ^. In the combinations of Eq. (fTUI) they 
are finite for z — > 1 and the imaginary parts stem exclusively from the complex logarithms. The 
arbitrariness in the representations of DVCS coefficient functions is exploited so that Eq. (j^Bjl re- 
duces to its equivalent in Eq. (jl3|) for 77 apart from the logarithmic divergences ~ log(m/j/ /i). 



5 Numerical results 

A suitable model for the GPDs to employ in a first numerical analysis of heavy quark effects is 
given by the moment-diagonal model defined in Ref. ^1]. One can expect that it leads to the same 
difficulties to describe experimental data, which is explained in detail in Ref. [IH], as factorized 
double distribution based models for GPDs [1^. Nevertheless moment diagonal models have 
important advantages. They are scale independently related to the conventional forward parton 
distribution functions (PDFs) so that no additional GPD input scale uncertainty is introduced 
and their dependence on the specific forward PDF set is reduced. Therefore, the following results 
should be unique consequences of heavy quark contributions. Apart from that, the evaluation 
of the form factors of the DVCS amplitude can be performed in a numerically very stable way 
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Figure 3: The imaginary and real part of the form factor Tj^, related to H in the appropriately 
substituted Eq. (j^, of the DVCS amplitude as a function of Q'^ at fixed ^ = 10"'^ and 
t = for three light flavors (solid line), with the fixed order charm contribution (dashed 
line), with the "massless" charm (dashed-dotted line) and additionally bottom (dotted line) 
contribution. 



without any principal value integrals |13]. The utilized unpolarized and polarized forward PDFs 
are GRV(98) [21] and the standard scenario of GRSV(OO) respectively, with three fixed light 
flavors. The charm and bottom quark masses are = 1.4 GeV and ruf, = 4.5 GeV. In accordance 
to (fT3j) and Ref. [21], intrinsic "massless" (off-)forward heavy quark distribution functions were 
generated from the same inputs. It should be noted that all the following fixed-order (FO) and 
"massless" parton (MP) results for the three light quarks u, d, s agree within 1.5%. 

In Fig. El the results for the form factor Tr^ ^{C, , t , Q"^) of the DVCS amplitude, which is the 
one that depends on H{x,^,t,fi), are presented as a function of for vanishing t and fixed 
^ = 10~^ which is representative for the kinematic region of DVCS measurements by HI [H] and 
by ZEUS [7_^ at DESY-HERA. The contributions of the intrinsic "massless" heavy quark GPDs 
show as expected a different threshold behavior compared to the fixed order results. The former 
start rather abrupt at = ml in the imaginary and real part whereas the latter set in smoothly 
at s > 4m^ in the imaginary part and contributes to the real part at any scale, though not 
sizeably for <C m\. Between the threshold region and the highest = 100 GeV^ displayed 
in Fig. El with the "massless" charm quark contribution in MP udsc lies significantly above the 
fixed order result FOudsc. The small "massless" bottom quark contribution justifies to neglect 
all further bottom results from now on, as it has been already done for the nearly vanishing fixed 
order bottom results in Fig. El 

The corresponding results with fixed C, = 10^^ relevant for the kinematic region of HERMES jl] 
and CLAS are shown in Fig. jH In this case the helicity dependent form factor Trp ^(^, t, Q^), 
which is the one that depends on H{x,^,t, fi) in the appropriately substituted Eq. Q, is also 
presented. At even smaller values of ^, e. g. ^ = 10"'^ in Fig. El Trp ^ becomes entirely negligible 
compared to Trj,^. The "massless" parton prescription still leads to considerable charm contri- 
butions whereas the fixed order charm results are marginal. Generally the corrections to DVCS 
observables due to heavy quarks get larger for increasing and in particular for decreasing ^. 
Therefore phenomenologically relevant effects have to be expected mainly for the DVCS cross 
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Figure 4: The imaginary and real parts of the form factors Tj ^ and TV a- related to H and H, 
respectively, of the DVCS amplitude as a function of at fixed ^ = 10~^ and i = for 
three light flavors (solid line), with the fixed order charm contribution (dashed line), and 
with the "massless" charm (dashed-dotted line) contribution. 
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Figure 5: The DVCS cross section as a function of Q"^ at fixed W = 89GeV for three light flavors 
(solid line), with the fixed order charm contribution (dashed line), and with the "massless" 
charm (dashed-dotted line), using a constant 6[GeV^^] = 4.9 as well as a strongly Q'^- 
dependent b{Q^)[Ge\/~^] = 4.9 - 5 log((5^[GeV]/100). The error bars denote the quadratic 
sums of the statistical and the systematic uncertainties of the ZEUS data 



section measurements at HI and ZEUS. 

The DVCS cross section was defined by using the equivalent photon approximation in Ref . . 
An approximate formula is given by jTHI 

aDvcs^^|TT,,(e,0,Q^)r (27) 

where 6 is a measure of the t-dependence of the DVCS amplitude and currently expected to 
be nearly constant. In Fig. the results for the DVCS cross section together with the ZEUS 
data at fixed W = \^ = 89 GeV are shown for a constant 6[GeV~^] = 4.9 as it has been 
used by the ZEUS collaboration and additionally for a highly Q"^ dependent 6((5^) [GeV~^] = 
4.9 — 5 log((5^[GeV]/100). The constant b shows the aforementioned difficulties by overshooting 
the data considerably that is basically in line with previous leading order predictions ^Uj- Only 
recently it was found that a "/c-delta ansatz" jlH] with k = 0, i. e. a narrow ?/-dependence 
for the double distribution F{x,y), can be used to describe the data However, such an 

ansatz is considered not realistic since it does not fulfill a certain symmetry constraint for double 
distributions [i^ liSj. 

The simple highly dependent 6-fit that yields a good agreement of the light flavor result 
with the data should merely give an impression of the relevance of charm contributions. The 
fixed order charm contributions increase the DVCS cross section by 10% for the lowest shown 
and up to 30% for large which is beyond the relative uncertainties of the first four data points. 
The "massless" charm GPD leads to significantly larger increases in the range of 60 . . . 90% which 
excessively overestimates the physical PCS subprocess. To check for the reliability of this result 
the CTEQ4 forward PDFs [20] have been used alternatively that give minor modifications less 
than 10% apart from a reduction of up to 15% around ~ 10 GeV^ that can be traced back to 
the higher charm mass tt^^teq = ^ g QeV and the correspondingly delayed generation of charm. 
Finally it should be noted that a change of the factorization scale fit a/ 4m^ + Q"^ changes the 
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fixed order results below the percent level. 

6 Conclusions and outlook 

The contributions of heavy quarks, especially of the charm quark, to deeply virtual Compton 
scattering have been adequately analyzed by fixed order perturbation theory via photon-gluon 
scattering. The use of intrinsic "massless" generalized heavy quark distributions overestimates 
the DVCS cross section at small skewedness < 10~^. In principle this may be absorbed in models 
for GPDs, which, however, is unreasonable since the fixed order calculation correctly describes at 
least the threshold region. At large scales both treatments cannot be distinguished because 
of present experimental uncertainties. 

The rough estimates indicate that heavy quark mass effects are relevant for small skewedness 
< 10~^ in particular at low scales with respect to experimental uncertainties. Nevertheless this 
should be confirmed with "reahstic" GPD models and using the full DVCS amplitude, possibly 
with higher twist corrections. For values of the skewedness ~ 0.1, characteristic for fixed-target 
experiments, contributions of heavy quarks can safely be neglected. The massive analytical results 
for the photon helicity fiip form factors of the DVCS amplitude have not been utilized but can 
be trivially included in future analyses. In this case, the massive description for heavy quarks 
should always be used as it smoothly reproduces a "massless" parton picture at high scales. 

To prove the pcrturbative stability of the PCS model the next-to-lowest order corrections to 
the massive hard scattering coefficients have to be calculated which are necessary for a consistent 
next-to-leading order analysis with heavy quarks correctly included. 
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